Document  Number:  SET  2015-0015 

412  TW-PA- 14245 


Estimators  for  iNET-Formatted  SOQPSK-TG 


May  2014 


Tom  Young 
SET  Executing  Agent 
412  TENG/ENI 
(661)  277-1071 

Email:  tommy.young.l@us.af.mil 


DISTRIBUTION  STATEMENT  A.  Approved  for  public  release:  distribution  unlimited. 


Test  Resource  Management  Center  (TRMC) 

Test  &  Evaluation/  Science  &  Technology  (T&E/S&T) 
Spectrum  Efficient  Technology  (SET) 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
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ABSTRACT 

This  paper  presents  algorithms  for  estimating  the  frequency  offset,  multipath  channel  coefficients, 
and  noise  variance  of  iNET-formatted  SOQPSK-TG.  The  estimators  compare  the  received  signal 
samples  corresponding  to  the  iNET  preamble  and  attached  sync  marker  (ASM)  bits  to  a  locally 
stored  copy  of  the  SOQPSK-TG  samples  corresponding  to  the  same.  The  mean  and  variance  of 
the  three  estimators  over  ten  test  channels  derived  from  channel  sounding  experiments  at  Edwards 
AFB  is  presented.  The  results  show  that  usable  estimates  are  achievable. 


INTRODUCTION 

Modem  digital  communication  techniques  leverage  the  capabilities  of  advanced  digital  signal  pro¬ 
cessing  capabilities  to  improve  the  BER  performance  of  modulated  carriers  in  harsh  radio  fre¬ 
quency  (RF)  environments.  A  fundamental  element  of  these  improvements  is  knowledge  of  the  RF 
environment  (e.g.,  frequency  offset,  multipath  channel  coefficients,  noise  variance).  If  the  trans¬ 
mitted  signal  contains  one  or  more  known  bit  fields,  then  a  sampled  data  receiver  may  compare 
the  received  signal  corresponding  to  the  known  bits  with  a  locally  stored  copy  of  the  known  bits  to 
estimate  the  desired  channel  conditions.  iNET-formatted  SOQPSK-TG  is  an  example  of  a  trans¬ 
mitted  signal  containing  a  known  bit  pattern  (in  the  form  of  preamble  and  ASM  bits:  see  the  top 
portion  of  Figure  1).  This  paper  describes  estimators  for  the  frequency  offset,  multipath  channel 
coefficients,  and  noise  variance. 

In  what  follows,  boldface  lowercase  variables  denote  column  vectors  whereas  boldface  uppercase 
variables  denote  matrices.  For  the  matrix  M,  MT  denotes  the  transpose  of  M  where  as  IVT  denotes 
the  Hermitian  (conjugate-transpose)  of  M. 
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The  transmitted  signal,  s(t),  is  SOQPSK-TG  whose  input  bit  stream  is  summarized  in  the  top 
portion  of  Figure  1.  The  signal  propagates  through  a  frequency  selective  channel  and  experiences 
a  frequency  offset  and  the  addition  of  additive  white  Gaussian  noise.  We  assume  the  received  signal 
is  filtered,  down-converted  to  complex  baseband,  and  sampled  (not  necessarily  in  that  order)  using 
standard  techniques.  The  resulting  sequence  of  received  samples  is 


r  ( n ) 


'  n2 

^  h(k)s{n  —  k ) 

_k=-Ni 


ej“on  +  w(n) 


(1) 


where  h(n)  is  the  impulse  response  of  the  unknown  channel  impulse  response  with  support  on 
— Ni  <  n  <  N2,  ojq  rads/sample  is  the  frequency  offset,  and  w{n)  is  a  proper  [1]  complex- valued 
zero-mean  Gaussian  random  process  with  auto-covariance  function 

{w{n)w*(n  -  k)}  =  <J2w8{k).  (2) 

As  an  example  of  the  need  for  these  estimators,  data-aided  MMSE  equalization  requires  estimates 
of  h(n)  for  —Ni  <  n  <  N2  and  [2].  But  good  estimates  for  h(n)  and  require  the  received 
data  to  be  de-rotated  by  an  estimate  of  the  frequency  offset  u0  rads/sample.  (A  good  estimate  for 
ujq  is  required  for  coherent  detection  as  well.)  The  estimator  described  in  this  paper  is  a  three  step 
estimator: 

1.  A  frequency  offset  estimate  u0  is  computed  from  the  received  samples  r(n)  corresponding 
to  the  iNET  preamble  bits. 

2.  The  channel  coefficient  estimates  h(n)  for  —N\  <  n  <  N2  are  computed  using  Cjq  and  the 
received  samples  r(n)  corresponding  to  the  iNET  preamble  and  ASM  bits. 

3.  The  noise  variance  estimate  is  computed  from  u>0,  h(n),  and  the  received  samples  r(n)  cor¬ 
responding  to  the  iNET  preamble  and  ASM  bits. 

Each  of  these  estimators  are  described  in  the  following  sections. 


FREQUENCY  OFFSET  ESTIMATOR 

Let  i  be  the  index  in  r(n)  corresponding  to  the  beginning  of  the  preamble  sequence  and  let  L,p 
be  the  number  of  samples  in  the  preamble  sequence.  The  situation  is  illustrated  in  Figure  1.  The 
complicating  factor  here  is  the  unknown  channel  impulse  response  h(n).  The  difficulty  of  the 
presence  of  h{n)  may  be  overcome  by  exploiting  the  periodic  properties  of  the  iNET  preamble. 

Because  the  iNet  preamble  comprises  8  repetitions  of  a  16  bit  sequence,  the  resulting  SOQPSK- 
TG  modulated  carrier  is  also  periodic  over  the  time  interval  corresponding  to  the  occurrence  of  the 
preamble.  Furthermore  the  sampled  SOQPSK-TG  signal  is  periodic  over  the  interval  correspond¬ 
ing  to  the  preamble  sequence.  The  situation  is  illustrated  in  Figure  1.  Here,  s(n)  are  the  samples 
of  the  complex  baseband  SOQPSK-TG  signal  and  p(n')  for  n!  —  0, 1, . . . ,  Lp  —  1  are  the  samples 
of  the  complex  baseband  SOQPSK-TG  signal  corresponding  to  the  preamble  bits.  If  the  preamble 


2 


CL  ^ 

2  ^ 

s  ^ 

.2  S 

a  o 
,3  c 

0)  <D 
Si  3 

0)  ^ 
rZ  CD 
33  c/d 

-S'  n 

S  o 

2  •* 

tz  cfl 

O  D 

■£>  -C 

^ '  -4— > 

«T  33 
33  3 

T3  ^ 

s  'e 


c/d 

< 

33 

3 

3 


?L 

c/d 

<D 


<D 


r-  ^ 


<D 


If 

H  & 
a  2 

«r  a 

■a  <u 

c  a 

"O  32 

<U  _ T 

■S  ^ 

O  cn 
«  <U 


2 

3 

•j— * 

O 

3 

b 


3 


3 
3 

<"  bX) 
•  *  • 

(U  cn 
M  a) 

2  -C 
ca  ■a 

&  _ 

IS 

33  L» 

H 
W 
Z  & 

•  s  £ 

<U  3 


<D 


33 

<D 

> 

’53 

o 

2a 


0) 


<D 

33  ^ 


Zh  ° 

3  C/D 
O  CD 
X 
<D 


3 

<D 

5-h 

CL 

<D 


2  c  w 

M 

^  <D 

<  .£2 


ID 


^  a 

<d  2 
2  2 
3  £ 

.SP  S 

Ph  J2 


a 

o 

o 

■4— > 

3 

33 


3 


sequence. 


sequence  starts  at  index  i,  then  s(i)  =  p( 0),  s(i  + 1)  =  p(  1)  and  so  on  to  s(i  +  LP  —  1)  =  p(Lp  —  1). 
Now  let  q(n)  be  the  samples  of  the  complex  baseband  SOQPSK-TG  signal  corresponding  to  the 
16-bit  sequence  that  is  repeated  to  form  the  preamble,  and  let  the  length  of  q{n)  be  Lq.  (At  an 
equivalent  sample  rate  of  2  samples/bit,  Lp  =  256  and  Lq  =  32.) 

Using  a  simple  change  of  index,  we  may  re-write  r(i),r(i  + 1), . . . ,  r(i  +  Lp  —  1)  as  r(i  +  £Lq  +  m) 
for  i  =  0, . . . ,  7  and  m  =  0. ....  Lq  —  1.  The  contribution  of  h(n)  to  r(i  +  tLq  +  m)  depends 
only  on  m,  the  position  inside  the  q(-)  block.  Consequently,  the  only  difference  (other  than  the 
additive  noise)  between  the  sample  at  position  m  in  adjacent  blocks  is  the  phase  rotation  due  to  the 
frequency  offset.1  This  feature  is  captured  by  the  correlation  function. 


Assuming  N\  <  Lq  and  N2  <  Lq,  we  are  interested  in  the  correlation  function  involving  the 
middle  six  blocks: 

i-\-7Lq  —  1 

R(S)  =  6L  _s  r(n)r*(n-S )  (3) 

^  n=i+Lq-\-8 

Because  the  samples  r(n)  are  periodic  with  period  Lq ,  we  are  interested  in  evaluating  this  function 
at  integer  multiples  of  Lq.  For  6  =  dLq  (1  <  d  <  5)  the  correlation  is 

2 

R(dLq)  =  ^-ejdL^°  +  v  (4) 

Lq 


where  a2  is  given  by 


Lq  —  1 


a2  =  ^  |a(m)|2 

m= 0 


(5) 


and  v  is  approximately  a  zero-mean  complex-valued  Gaussian  random  variable.  Inspired  by  the 
delay-and-multiply  (D&M)  estimator  described  by  Eq.  (4.5.1)  of  [3],  we  use  a  frequency  estimator 
based  on  (4)  corresponding  to  d  —  1: 


(6) 


This  estimator  is  simple  and  exploits  the  noise  averaging  that  occurs  over  the  5 Lq  consecutive 
samples  of  r(n). 


CHANNEL  ESTIMATOR 

The  channel  estimator  is  based  on  the  SOQPSK-TG  samples  corresponding  to  the  iNET  preamble 
and  ASM  fields.2  The  situation  illustrated  in  Figure  2.  Because  the  channel  estimator  does  not 

'This  is  only  true  for  middle  six  blocks  (i.e.,  £  =  1 .....  6)  because  the  convolution  sum  in  the  first  block  (£  =  0) 
includes  samples  from  the  data  field  immediately  preceding  the  first  block  whereas  the  last  block  (£  =  7)  includes 
samples  from  the  ASM  field  immediately  following  the  last  block.  Eliminating  the  first  and  last  blocks  is  sufficient  as 
long  as  Ni  <  Lq  and  N2  <  Lq.  If  either  of  these  conditions  are  not  true,  then  this  is  true  for  fewer  blocks. 

2Because  the  preamble  is  an  8-times  repetition  of  CD98|lex.  the  convolution  matrix  formed  from  the  samples  cor¬ 
responding  to  the  preamble  is  not  sufficiently  well-conditioned  to  provide  robust  channel  estimates  in  the  presence  of 


4 
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exploit  the  periodic  nature  of  the  preamble,  the  periodic  elements  of  the  preamble  sequence  are  not 
made  explicit  here.  For  notational  purposes,  we  let  x(0), . . . ,  x(Lp  +  La  —  1)  represent  the  samples 
corresponding  to  the  preamble  followed  by  the  ASM  bits. 

The  first  step  is  to  remove  the  frequency  offset  from  the  sequence  of  received  samples.  This  is 
accomplished  using 

rd(n)  =  r(n)e~JU}°n,  n  =  i, . . .  ,i  +  Lp  +  La  +  Ld  -  1  (7) 

where  Ld  =  N  x  6144  is  the  number  of  samples  corresponding  to  the  6144  data  bits.  (Here  the 
subscript  d  means  “de-rotation.”) 

The  next  step  is  to  organize  the  rd(n )  corresponding  to  the  preamble  and  ASM  bits  in  to  a  vector, 
but  skipping  the  first  N2  samples  and  last  Nl  samples: 

rd(i  +  N2) 

rd=  !  .  (8) 

r  d{i  +  Lp  +  La  —  Ni  —  1) 

This  vector  is  related  to  the  samples  corresponding  to  the  preamble  and  ASM  samples  by 

rd  =  Xh  +  w  (9) 

where  X  be  the  (Np  +  Na  —  Ni  —  N2)  x  ( Ni  +  N2  +  1)  convolution  matrix  formed  from  these 
samples: 

x(N2  +  Ni)  ■  ■  ■  z(0) 

X  =  :  :  (10) 

x(Lp  +  La  —  1)  •  •  •  x(Lp  +  La  —  Ni  —  N 2  —  1) 

and  where  h  is  the  ( Ni  +  +  1)  x  1  vector  formed  by  the  channel  coefficients, 

h  =[h(-Ni)  •••  h(N2)) T  ,  (11) 

and  w  is  the  (Lp  +  La  —  Ar,  —  N2)  x  1  vector  of  noise  samples.  From  this,  the  vector  of  channel 
estimates  is 

h=  (Xtx)-1Xtrd.  (12) 

At  first  glance,  the  computational  complexity  of  (12)  seems  quite  high.  But  note  that  the  matrix 
X  is  fixed.  Consequently,  the  matrix  product  (X^X)  1  X^  may  be  precomputed  and  stored  at  the 
receiver.  The  result  is  that  the  channel  estimator  reduces  to  a  matrix- vector  multiplication. 

even  a  small  frequency  offset.  The  convolution  matrix  formed  from  samples  corresponding  to  both  the  preamble  and 
ASM  fields  has  less  structure  and  provides  more  robust  channel  estimates.  See  Section  4  of  [4], 
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NOISE  VARIANCE  ESTIMATOR 


The  noise  variance  may  be  computed  from  a  scaled  version  of  the  sample  variance  corresponding  to 
the  error  between  the  observed  signal  and  its  estimate.  In  matrix/vector  form,  the  sample  variance 
is  the  inner  product  of  a  vector,  divided  by  its  length.  The  noise  variance  estimator  is 


a?„  = 


2  P 


rd-Xh 


(13) 


where  rd  is  given  by  (9),  X  by  (10),  and  h  by  (12).  The  scale  factor  p  is 

p  =  Tr{l-X(XtX)“1Xt}  (14) 

where  I  is  the  (Np  +  Na  —  Ni  —  N2)  x  (Np  +  Na  —  Ni  —  N2)  identity  matrix.  This  scale  factor  is 
chosen  to  produce  an  unbiased  estimate  of  the  noise  variance  (see  Section  7  of  [4]). 


PERFORMANCE 

The  performance  of  the  three  estimators  was  evaluated  in  compute  simulation  using  the  following 
parameters: 

1 .  The  payload  data  rate  was  equivalent  to  10  Mbits/s  (the  equivalent  “over-the-air”  bit  rate  was 
10.3125  Mbits/s).  The  iNET-formatted  SOQPSK-TG  signal  and  channel  were  generated  at 
an  equivalent  sample  rate  of  2  samples/bit. 

2.  Because  the  channel  estimator  does  not  know  the  true  length  of  the  channel,  the  estimator 
used  values  for  IV i  and  N2  larger  than  any  of  the  test  channels.  These  values  were  Ni  =  12 
and  N2  =  25  samples. 

3 .  The  simulations  were  performed  over  10  representative  channels  derived  from  channel  sound¬ 
ing  measurements  conducted  at  Edwards  AFB  under  the  M4A  program  [5],  The  test  chan¬ 
nels  are  summarized  in  Table  1  and  the  corresponding  frequency-domain  plots  are  shown  in 
Figure  3. 

The  results  for  the  frequency  offset  estimator  (6)  are  shown  in  Figures  4  and  5.  The  simulated 
estimator  means  for  frequency  offsets  —50  kHz  <  A/  <  50  kHz  over  the  ten  test  channels  are 
plotted  in  Figure  4.  This  figure  shows  that  the  frequency  estimator  is  unbiased  and  capable  of 
accurately  estimating  the  frequency  offset  over  this  range.  The  simulated  variance  of  the  frequency 
offset  estimator  (6)  is  shown  in  Figure  5  where  it  is  observed  that  the  estimator  error  variance 
decreases  with  increasing  Eb/N0,  as  expected.  The  RMS  frequency  error  is  the  square  root  of  the 
variance.  The  plots  show  that  for  Eb/N0  >  12  dB,  the  RMS  frequency  error  is  a  few  hundred  Hz. 
This  is  small,  but  a  mechanism  for  tracking  out  the  residual  frequency  offset  in  the  SOQPSK-TG 
detector  is  still  required.  (For  example,  see  Figure  3  of  [2]). 
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Table  1:  Description  of  the  ten  test  channels  used  in  the  simulations. 


channel 

iV, 

n2 

length 

environment 

01 

2  1 

7 

9 

Taxiway  E 

02 

2 

17 

20 

Taxiway  E 

03 

1 

22 

24 

Taxiway  E 

04 

5 

13 

19 

Takeoff  on  22L 

05 

1 

1 

3 

Cords  Road 

06 

1 

2 

4 

Cords  Road 

07 

0 

4 

5 

Cords  Road 

08 

2 

3 

6 

Black  Mountain 

09 

1 

1 

3 

Black  Mountain 

10 

2 

3 

6 

Land  on  22L 

The  results  for  the  channel  estimator  (12)  are  shown  in  Figures  6  and  7.  Here,  the  figure  of  merit 
is  based  on  the  sum  of  estimation  errors: 


n2 

£h=  [ h(n )  -  h{n) 

n=—N\ 


(15) 


In  Figure  6,  the  error  mean  is  the  expected  value  of  £/-,.  which,  in  turn,  was  estimated  in  simulation 
using  the  sample  mean.  In  Figure  7,  the  error  variance  is  the  variance  of  £/,  which,  in  turn,  was 
estimated  in  simulation  using  the  sample  variance.  The  impact  of  the  frequency  offset  estimation 
errors  on  the  accuracy  of  the  channel  estimates  is  clearly  evident.  The  estimator  error  mean  only 
decreases  with  increasing  Eb/N0  when  Eb/N0  >  15  dB. 


The  results  for  the  noise  variance  estimator  (13)  are  shown  in  Figures  8  and  9.  Because  the  noise 
variance  decreases  with  increasing  Eb/N0,  we  plot  the  normalized  versions  of  the  noise  variance 
estimator  error: 


£  = 


02 


<72 


of, 


(16) 


As  with  the  channel  estimator,  the  noise  variance  estimator  clearly  suffers  from  the  frequency  and 
channel  estimation  errors.  Consequently,  of  the  three  estimates,  this  one  is  least  reliable. 


How  these  estimation  errors  impact  the  overall  bit  error  rate  (BER)  is  discussed  in  detail  in  [4] .  The 
impact  of  the  estimation  errors  may  also  be  inferred  from  the  (BER)  performance  of  the  MMSE 
equalizer  presented  in  [2]  where  it  is  demonstrated  that  the  equalization  techniques  enabled  by 
these  estimates  have  reasonably  good  BER  performance  over  the  same  set  of  ten  test  channels. 


CONCLUSIONS 

This  paper  presents  algorithms  for  estimating  the  frequency  offset  (6),  multipath  channel  coeffi¬ 
cients  (12),  and  noise  variance  (13)  of  iNET-formatted  SOQPSK-TG.  The  estimators  compare  the 
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received  signal  samples  corresponding  to  the  iNET  preamble  and  attached  sync  marker  (ASM)  bits 
to  a  locally  stored  copy  of  the  SOQPSK-TG  samples  corresponding  to  the  same.  The  mean  and 
variance  of  the  three  estimators  over  ten  test  channels  derived  from  channel  sounding  experiments 
at  Edwards  AFB  is  presented.  The  impact  of  the  estimation  errors  may  be  inferred  from  the  (BER) 
performance  of  the  MMSE  equalizer  presented  in  [2]  where  it  is  demonstrated  that  the  equalization 
techniques  enabled  by  these  estimates  have  reasonably  good  BER  performance  over  the  same  set 
of  ten  test  channels. 
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Figure  3:  Frequency-domain  plots  of  the  example  channels  from  channel  sounding  experiments  at 
Edwards  AFB.  In  each  plot,  the  thick  line  is  the  channel  frequency  response  and  the  thin  line  is  the 
power  spectral  density  of  SOQPSK-TG  operating  at  10.3125  Mbits/s. 


10 


estimator  mean  estimator  mean  estimator  mean  estimator  mean  estimator  mean 


-50* 


0  50 

true  frequency  (kHz) 


-50* 


0  50 

true  frequency  (kHz) 


-50* 


0  50 

true  frequency  (kHz) 


-50* 


0  50 

true  frequency  (kHz) 


-50* 

-50 


0  50 

true  frequency  (kHz) 


to 

a> 

E 

o 

to 

E 


-50* 


-50  0  50 

true  frequency  (kHz) 


03 

0 

E 

b 

C0 

E 


-50  0  50 

true  frequency  (kHz) 


03 

0 

E 

b 

0 

E 


-50  0  50 

true  frequency  (kHz) 


CO 

a> 

E 

o 

to 

E 


-50  0  50 

true  frequency  (kHz) 


50 


CO 

CD 

E 

o 

to 

E 


-50 


test  channel  10 

- ^ 

jm 

m" 

,« 

-* 

Jtr 

-50  0  50 

true  frequency  (kHz) 


Figure  4:  Estimator  means  for  the  frequency  estimator  (6)  for  the  ten  test  channels  described  in 
Figure  3. 
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Figure  5:  Estimator  error  variance  for  the  frequency  estimator  (6)  for  the  ten  test  channels  described 
in  Figure  3. 
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Figure  6:  Estimator  means  for  the  channel  estimator  (12)  for  the  ten  test  channels  described  in 
Figure  3. 
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Figure  7:  Estimator  error  variance  for  the  channel  estimator  (12)  for  the  ten  test  channels  described 
in  Figure  3. 
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Figure  8:  Estimator  means  for  the  noise  variance  estimator  (13)  for  the  ten  test  channels  described 
in  Figure  3. 
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Figure  9:  Estimator  error  variance  for  the  noise  variance  estimator  (13)  for  the  ten  test  channels 
described  in  Figure  3. 
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